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Abstract 



A predictive control scheme for a permanent-magnet synchronous machine (PMSM) is presented. It is 
based on a suboptimal method for computationally efficient trajectory generation based on continuous 
kj '. parameterization and linear programming. The torque controller optimizes a quadratic cost consist- 

■ mg of control error and machine losses in real-time respecting voltage and cuiTcnt limitations. The 

^ I multivariable controller decouples the two cuiTcnt components and exploits cross-coupling effects in the 

c/3 ! long-range constrained predictive control strategy. The optimization results in fast and smooth torque dy- 

CJ ■ namics while inherently using field-weakening to improve the power efficiency and the current dynamics 

in high speed operation. The performance of the scheme is demonstrated by experimental results. 

^ ■ Introduction 

O ■ The efforts of implementing predictive controllers in electrical drives aim at replacing the classical 

O i cascaded field-oriented control structure with PI controllers. The machine can be better exploited by 

_i ! improved control behavior, the system variables are optimized. In this contribution, the conventional 

^ torque and current control structure of two separate controllers is changed to multi-input multi-output 

(MIMO) control. By transformation into the field-oriented frame, torque generation is decoupled from 
flux variation, however, the current controllers are still strongly coupled, therefore a MIMO controller is 
advantageous. The obtained improvements are better decoupling, better cuiTcnt and voltage constraint 
handling by exploiting cross-coupling between the orthogonal components, and better power efficiency 
and dynamics by optimally adjusting the currents in both dynamic and steady-state operation. 
The major obstacle in implementing predictive control schemes is the limited computational power, in- 
i herited by the high sampling rates. Existing implementations suffer from this restriction and can not 

exploit the full advantages of model predictive control (MPC). For instance, generalized predictive con- 
trol (GPC) has a high optimization horizon but is unconstrained, whereas predictive torque control (PTC) 
is constrained but so far only reaches 2 steps of prediction [T|. For good performance, both, inclusion of 
constraints and a high optimization horizon are required. Using control in the field-oriented frame, the 
analytical problem description enables using mathematical optimization algorithms. Such schemes are 
computationally efficient and maximize the obtained information for a given computational power 
The online solution of the linearly constrained linear-quadratic problem, typical for MPC, requires 
quadratic programming (QP) algorithms, which are, however, computationally too expensive for drive 
systems. A recent development is the use of explicit MPC, where an offline solution is computed and 
stored as look-up table in the real-time controller |]2l Q. The scheme reaches 5 prediction steps with 
constraints. While several fast online algorithms were recently proposed |4| [5 |, most of them based on 
gradient search algorithms, the interest in online solutions is growing. Advantages of online optimization 
are manipulation and adaptation of parameters, resulting in a more flexible choice of machines for the 
control system. 

This implementation is based on a suboptimal trajectory generation algorithm presented in [61, embedded 
in a flatness-based predictive control scheme Q- It is based on a continuous approach, the variables 
are not discretized but represented as a polynomial with undetermined coefficients. With this method, 
higher optimization horizons can be reached with comparably few parameters, considerably reducing 
the computational burden. Constraints are handled by linearization of the cost functional and the use 
of a lineal- programming (LP) solver, which is amongst the smallest and fastest numerical optimizers. 



As result, a (suboptimal) prediction of 2 ms with current and voltage constraints is obtained at 8 kHz 
sampling rate. 



Problem Statement 

Machine Model 

As linear-quadratic optimization problems with linear constraints are simpler to solve in real-time, the 
machine model is linearized. Assuming that the rotor speed does not change too much over the optimiza- 
tion horizon T , 

—(OM(t)^O ^ (Om(0 = const. V?G[0,rl, (1) 
at 

the PMSM model and the voltage equations become linear. The electrical subsystem of the machine, 
consisting of the quadrature and direct currents and id (peak values), is given as 

d 

^d-^^id = -Rid + npdiMLqiq + ^d, (2) 



'^^dt'^ = -Riq - npiOnfLdid - fipiOMK + Uq, (3) 
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'Cm = 7;npKiq. (4) 
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The nomenclature is shown in table II in the appendix. It is noted that the reluctance torque = 
|?ip(Lj —Lq)iiiiq is neglected, as this term is very small compai^ed to the electromagnetic torque in 
surface-mounted PMSMs or in machines with small saliency. Furthermore, it would render the model 
nonlinear, requiring nonlinear optimization methods [Si. 



Optimization Goals and Cost Functional 

The formulation of a suitable cost functional is a key point in predictive control, as it is the only tuning 
possibility of the control scheme. The optimization is aiming at minimizing the control enor for good 
dynamical performance as well as machme losses for better efficiency. Both goals are included in the 
cost functional. By choosing the cost functional and weights well, it is possible to find a good trade- 
off between both goals during transients, or eventually to fulfill both goals in steady-state. The cost 
functional for the predictive torque controller is 

J = C {Pctrl{t) + Wl ■ Ploss{t)) dt + T- P,,rl{T) (5) 

Jo 

which trades off the squared control eiTor 

Pctrl{t) = {^M-^*M)^ (6) 

with machine losses 

PioAt) = ^Riid + il) + \f'p^MkFei'i'd + 'i'l) = ^R{il + il) + ^np(OMkFe{iLdid+Kf + iLqiq)^). 

(V) 

The first term in Pjoss represents copper losses, and the second term represents the iron losses consisting 
of hysteresis losses. Eddy current losses are negligible on the tested machines, however, they could 
be included using the model presented in [9]. Iron losses can be reduced by field- weakening, where a 
trade-off between copper and iron losses is found lITOl . Imposing a negative direct cuiTcnt i^, the flux 
magnitude in the stator is reduced while the copper losses increase. As the iron loss constant kfe is not 
part of standard motor parameters, it has to be determined experimentally ifTTTl . The last term in / is the 
end-weight of the control eiTor and aims at reducing the steady-state control error. The weight Wl was 
set 0.05, the value was determined heuristically. 

The optimization horizon is set T = 2 ms such that the cost functional includes the complete setpoint 
change. It is important that the optimization horizon is high enough, otherwise the open-loop and closed- 
loop trajectories differ and the behavior is strongly suboptimal. This is illustrated in Fig. [T] If the 
horizon is too small, due to the end-weight of the control error, a significant difference between the 
optimized open-loop trajectories, and the closed-loop trajectories resulting from regeneration at every 
sampling step, appears. Then, the closed-loop trajectories simply don't fit the cost functional anymore. 
For a horizon higher than required for the setpoint change, the difference between open- and closed-loop 
trajectories becomes smaller. 



Figure 1: Open- and closed-loop trajectories in predictive control. Left: small horizon, right: high horizon. 



Current and Voltage Constraints 

The most important nonlinearities of a PMSM, in terms of control, are the voltage and current limitations. 
The cuiTcnt constraints prevent overheating of the machine, and the voltage is limited by the maximum 
DC-link voltage of the inverter. The voltage constraints limit rotor speed as well as current dynamics in 
high-speed operation. These constraints are linearized, in order to be computationally efficiently treated. 
The current range for the direct current id is limited to ;™"' < < 0. Only negative values of are desir- 
able, as they improve power efficiency and reduce the induced voltage by weakening the flux magnitude 
in the stator [31 ifTOl . The lowest value i""" is the optimum value at rated speed {-^Pioss = 0) and is given 
as 
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2 I R 



(8) 



which is independent of quatrature current iq as the reluctance torque was neglected. The value is dou- 
bled to enable further field-weakening to improve dynamics in high speed, an effect described in the 
experimental results section. For the quadrature current iq, the largest possible range of values should be 
available. The resulting linear constraints, shown on Fig. |2j almost completely fill the current region of 
interest. A linearization is thus acceptable. 

The voltage linearization is a little bit more difficult. The ^-axis should not be restricted, as the induced 
voltage is aligned to it and is the largest value that will appear. A steady-state analysis of the system 
equations (2), (3) shows that a rectangular voltage area results 

RiT - npLqiori'r ^ ^ npLqCorir, (9) 
-Ri'j;"' + npLd(0Tit - npKdr < Uq < Ri'q""' + npK(oZ'\ (10) 

This rectangle is expanded such that the outer circle of the voltage limitation is hit (light grey on Fig. 
During dynamical transients, the voltage vector points to one of the outer corners, subsequently touching 
the outer limiting circle. Therefore, linearizing the voltage limits as a rectangle by the presented method, 
as shown on Fig. [2l does not limit the operational range and only marginally affects dynamics. A 
less restrictive method is presented in 0, where a time- varying linearization in form of a hexagon in 
stator frame is presented. While such a linearization is possible with the underlying predictive control 
algorithm, the chosen method in the (<i,^)-frame is chosen for simplicity. 



Optimal Control Algorithm 

To study real-time applicability of the presented scheme, first, the highest possible amount of optimiza- 
tion parameters is determined. The fastest optimizer with constraints is the widely known linear pro- 
gramming (LP) method. Table I shows some worst-case computational results of LP (simplex method 
from [|12|]) as function of the number of free parameters (CPU: 1.4 GHz industrial PC). More parameters 
lead to a higher number of iterations which are also more complex; the worst-case number of iterations 
is the number of parameters plus the number of constraints. As in the underlying application, the con- 
straints are decoupled, however, this worst-case is not to be expected. The maximum runtime is given 
by the sampling rate minus latency of input/output, therefore at 8 kHz sampling rate, it must be less than 
about 110 /us. Thus, at best, 12 parameters can be optimized. Runtime of the predictive controller is 
further discussed at the end of the section. 



Figure 2: Linearized current and voltage constraints. Circle: feasible set of current and voltage vectors, grey: 
feasible set after linearization of the constraints. 



Table 1: Runtime of a linear program for some worst-case problems on a 1 .4 GHz CPU 



Parameters 


Constraints 


Iterations 


Runtime [ps] 


20 


48 


67 


769 


12 


32 


34 


165 


8 


14 


10 


35 



Trajectory Generation 

The trajectory generation algorithm presented in f6l, a development related to flatness-based meth- 
ods |[T3l . was designed for this application. It can optimize a quadratical cost function with linear 
constraints. As major differences to standard algorithms, it is applying a continuous parameterization 
instead discretization, and the computationally efficient linear programming solver is used instead of 
quadratic programming or iterative gradient search. 

The trajectories for the cuiTcnt are defined as degree n power series with undetermined coefficients a,y, 

n n j.k 

id{t) = £ adk ^, iq{t) =Y.^qk7fJl, te [0,T]. (11) 
k=0 ^ k=0 ^ 

Due to the analyzed computational limitations, n = 3 is chosen as polynomial degree. The first co- 
efficients a^o and a^yo are the initial conditions, and the remaining 6 coefficients are determined by 
optimization. The corresponding voltages Uja{t) are computed by algebraic differentiation of (fTTI) and 
by solving for the model equations Q and ©, this is also called the flatness-based approach |[T4l . This 
way, the voltages do not need to be represented by additional parameters. 

Substituting the variables i^q, ^^idq and Udq by the found functionals in the cost functional J Q and 

constraints can easily be done with a computer-algebra tool (for instance Maplesoft's® Maple^'^). The 
cost functional / is then a quadratic function of the unknown parameters, of the machine parameters, the 
measured currents, the speed ©m and the torque reference x^. Graphically, J can be represented as in 
Fig. [3] (left). 

As / is convex, the unconstrained optimum is found algebraically by solving first-order necessary con- 
ditions, in this case, the solution of a system of linear equations by matrix inversion. Then, by an affine 
coordinate transformation, the problem can be reformulated as least-distance problem, i.e. a quadratical 
cost describing the distance to the unconstrained optimum. As result, the cost functional looks much 
simpler, see Fig. [3] (middle). 

In the next step, the least-distance problem is linearized around the unconstrained optimum, see Fig. [3] 
(right). As the coordinate transformation to the least-distance problem is an affine transformation, and as 
the linearization doesn't affect the constraints, it is obvious that the linear constraints remain linear after 
this transformation. One point is, however, more difficult, the parameterization of the constraints with the 
polynomial coefficients, presented in [6]. The linearization of the cost function inherits a large eiTor in 
the value of /, but the values of the coefficients a are not affected that much: the least-distance problem 
is not so much different in the linear form as it would have been in quadratical form. Furthermore, see 
that a difference only appears if a constraint is active, the unconstrained optimum is the same. It can be 
shown that the resulting cost inherited by the linearization is 

j' = jQ + 2n-Jc (12) 

in the worst case, where Jq is the unconstrained cost, and Jc the extra cost when considering constraints. 
The suboptimality is therefore bounded. As « = 3, there are 6 parameters in the quadratical cost function. 




Figure 3: Trajectory generation algorithm. Left: original problem, middle: transformed least-distance problem, 
right: transformed and linearized problem. 



During the linearization, this number is increased to 12, as in the LP standard form, only positive values 
are possible. After the transformation, the problem is reformulated in standard form for linear program- 
ming, and a simplex solver [12:] can be run. To visualize the suboptimality, a comparison (simulations) 
between the linearized and the original problem is presented in [6J. 

All presented computations are done using Maplesoft's® Maple^'^, from where the matrices for the 
LP solver are generated using a C code generation toolbox. The real-time software thus consists of 
a simplex tableau assignment, which is automatically generated code, a simplex LP solver and some 
post-processing. As the assignment is based on symbolic calculations, the machine parameters can be 
changed online. 



Predictive Control 

The trajectory generation scheme is embedded in a predictive controller. The control structure is shown 
in Fig. |4l A cascaded control structure is chosen as speed is assumed constant for trajectory generation. 
As the mechanical plant is generally only roughly known, this structure is advantageous. Model-based 
control can be used for the machine as the parameters are known, but for the mechanical part, any robust 
feedback controller can be chosen. 
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Figure 4: Control structure of the predictive torque controller cascaded by PI speed control. 



First, measurement and control timing is analyzed. Two timing sequences are shown on Fig. [5] (a) and 
(b). The interrupt-based control system triggers an interrupt every 125/js (green signal). At this instant, 
the applied voltage command (magenta signal) is modulated by space-vector modulation. At the same 
instant, the A/D conversion of the cuiTcnt measurements (blue signal) is performed to avoid the impact of 
current ripples. After the interrupt, the current iq [k] is available and the controller (cyan signal) is started 
to compute the next voltage command Uq[k + \], to be modulated at the next interrupt. The computational 
delay is accounted with a delay compensation technique that is applied to generate iq[k + 1] as feedback 
value for the control law to compute + 1] lITSl . it corresponds to a prediction of one sampling step 
with the model, the cuiTcnt iq[k] as well as the previously commanded voltage u[k]. Furthermore, see that 



the response to the commanded voltage Uq[k] is iq[k+ 1] which is available one interrupt later, this is the 
plant delay and it is naturally included by recalculating the trajectory at every sampling step, see Fig. [T] 
From the predicted trajectory, the applied voltage is Udq[k] = Udq{0). 

From the timing sequences, interesting insight into the computational demands of the algorithm is gained. 
The first part of the cyan signal shows the calculation time for the simplex tableau initialization, it takes 
about 10/JS. Included in these calculations, which are results of symbolical calculations, is a calculation 
of the unconstrained optimum and the linearization of the problem. The second and biggest part of the 
cyan signal is the runtime of the linear program. At the beginning of Fig. [5] (a), where voltage and 
current are both zero, it is only about 20/is, but to calculate the voltage step at 2000 rpm shown on Fig. 
|5] (b), more iterations are involved as many constraints are active, and the computation time rises to 
almost 60/JS. The total time of the inteiTupt handling, the simplex initialization, the LP solver and the 
post-processing sum up to almost 100;us in the worst case, therefore up to 80% of the available time is 
used. 



Results 

A surface-mounted PMSM with parameters shown in table II is used. It is coupled to a load drive such 
that an arbitrary load torque or speed can be applied. The algorithm is implemented on a PC- 104 based 
real-time system with a 1.4 GHz CPU described in |[T6ll . The voltage limitation was set to 75% of the 
possible 330 V to clearly see behavior of saturation in control, accordingly, rated speed is reduced from 
3000 to 2250 rpm. 

Experimental results of the proposed scheme are shown in Fig. [51 Subfigure (c) shows the response 
to two subsequent speed reference steps, the load drive is deactivated. The torque is increased rapidly 
by a high but feasible voltage peak. Interestingly, no overshoot arises, even though the speed change is 
performed quickly. The direct current is proportional to the speed and thereby reduces iron losses which 
are considerable at high speeds. Losses are decreased by about 4%, and the efficiency is improved by 
about 0.5% at 2000 rpm. Better results are obtained on machines with higher inductances ifTOll . 
The next three subfigures (d), (e) and (f) show fast torque transients at zero, medium and high speed, 
respectively. The PMSM is in torque control mode while the load drive keeps speed constant at 0, 2000 
and 2400 rpm, respectively. The current components are well decoupled, a fast current change on the 
quadrature axis does not affect the direct axis at all in (d) and (e). With two separate PI controllers, a 
short current excursion would be seen on the direct axis during the torque transient. Again, the current 
on the direct axis /V/ is proportional to the speed. Furthermore, the torque is fast and at the same time 
smooth, the voltage becomes proportionally small for smaller control errors - a nice characteristic of 
quadratic cost functional, compared to linear cost functions which result in deadbeat behavior and are 
more sensitive to uncertainties llT5l . 

On subfigures (d) and (e), behavior with active voltage constraint is the same as when using standard 
saturation or anti-windup strategies. On subfigure (f), however, a different behavior is seen, the direct 
current id is reduced to perform field-weakening. This implies that the stator induced voltage is reduced 

on the quadrature axis, see eq. Q. Thereby the derivative of quadrature current ^Jq is higher and the 
torque-generation dynamics increased, at the cost of higher copper losses on the direct axis. Without 
additional field-weakening, the reference torque would not be reached after the optimization horizon 
of 2 ms, thereby the end-weight in / oversizes the loss term. Therefore, in this predictive control im- 
plementation, field-weakening not only improves efficiency, but also improves dynamics by exploiting 
cross-coupling between the orthogonal cuiTcnt components to optimally bypass the voltage saturation. 
It is also possible to operate the PMSM beyond rated speed by steadily doing field-weakening to bypass 
the voltage saturation on the quadrature axis, as shown in [31- It is remarked that the current on the direct 
axis id has no reference, its value is obtained from the optimization of the cost function. Therefore the 
method works well and is numerically stable; the optimal value follows inherently. 



Conclusion 

A predictive control scheme for a PMSM was introduced. Based on suboptimal real-time optimization, 
the currents and voltages are computed according to a cost functional. The prediction horizon is 2 ms 
at a sampling rate of 8 kHz, and voltage and cun^ent constraints ai^e respected. The advantages of the 
long-range constrained predictive MIMO control scheme can be concluded as follows: improved decou- 
pling and accounting of cross-coupling, precise measurement and control timing, respecting current and 
voltage constraints, fast and smooth dynamical behavior, improved power efficiency by field weakening, 
and improved dynamics close to voltage saturation by field-weakenmg. 

Furthermore, it was shown that it is possible to implement long-range MPC with constrained online- 
optimization even on fast-sampling systems such as electrical drives. 
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Fi gure 5: Experimental results of the predictive control scheme. Subfigures (a)-(b): blue: quadrature current iq 
without delay compensation (5.3 A/div), cyan: quadrature voltage (160 V/div), green: interrupt handling, cyan: 
control law computation. Subfigures (c)-(f): green: rotor speed com (1500 rpm/div), blue: quadrature current iq 
with compensation (5.3 A/div), cyan: quadrature voltage Uq (160 V/div), cyan: direct current (2.1 A/div). 



Appendix: Machine Parameters 



Table II: Nominal Parameters of the Synchronous Machine 



IVTnTiiifnr'tiirpr ^ IVToHpI 


Merkes MT5 1 050 


Rated Power Pf^ 


2640 W 


Rated Torque Xmn 


8.4 Nm 


Rated CuiTent (peak) 


5.6 A 


Rated Speed (Hmn 


3000 ipm 


Pole Pairs Hp 


3 


Rated Voltage (peak) 


560 V 


Stator Inductance L^, Lq 


4.8, 7.2 mH 


Stator Resistance 7? 


0.92 n 


Motor Constant K (peak) 


0.334 Vs 


Iron Loss Constant (Hysteresis) kpe 


1.27 A 

Vs 
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